Extending holographic LEED to ordered small unit cell superstructures 
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I. INTRODUCTION 

The short inelastic scattering length of low energy elec- 
trons makes them ideal tools for surface crystallography. 
The energy variations of the intensities (I(V) curves) 
of the reflected Bragg beams depend sensitively on the 
structure of the few outermost atomic layers, which is 
usually deduced by a trial-and-error fitting to simulated 
I(V) curves from model structures^. 

Even for a single model structure the calculation of 
a set of I(V) curves is no trivial task. The calculation 
requires two main ingredients: first, a representation of 
the scattering properties of each of the atoms of the pre- 
sumed structure - this is parametrized in a set of energy- 
dependent phase shifts for each angular momentum quan- 
tum number. Second, it is necessary to evaluate the scat- 
tering of electrons amongst all atoms of the sample. Dur- 
ing the 1960's and 1970's many ingenious methods were 
developed to efficiently sum the multiple-scatterirujupaths 
followed by the electrons from source to detectorEM 

Although many of the simpler surface structures were 
solved during this period, it was also becoming clear that 
more complicated structures would probably forever lie 
beyond the reach of standard multiple-scattering trial- 
and-error fitting methods due to the exponential scaling 
of computer time with the number of parameters to be fit- 
ted. A simple example makes this abundantly clear: sup- 
pose we need to determine the 3 Cartesian coordinates of 
N symmetry inequivalent atoms, and for this purpose we 
consider models with 10 values of each coordinate. The 
time required for an exhaustive search amongst these pa- 
rameters is T = S x 10 3N , where S is the time required 
to calculate the LEED spectra from each model. Taking 
S = 1 second and N = 3 we find T = 30 years. Just 
doubling the number of atoms to N — 6 results in a time 
T of the order of the age of the universe. In the latter 
case, increases in computer speeds by even a million-fold 



would only bring T down to a few millenia! (Of course, 
the above time estimates reduce considerably if it were 
judged that not all three Cartesian coordinates need to 
be determined independently). 

Obviously there is a need to examine some alternatives 
to exhaustive trial-and-error searches. Consequently, 
procedures applying a directed search in the multidimen- 
sional parameter space have been developed (for reviews 
see e.g refs. ||,|]). In this context the testing of simulated 
annealing-like and genetic algorithmsU are certainly steps 
in a most promising direction. Also, quasidirect meth- 
ods were developed allowing a search-free optimization 
of structural parameters once a model near the correct 
structure had been found or guessedB. The other ap- 
proach that has shown considerable promise for_LEED 
in recent years has been the holographic methodEI. 

This has been applied primarily to surfaces containing 
adatoms of either the same or different species than the 
surface atoms. Initially it was used to reconstruct a 3D 
image of the local geometry of disordered atomic adsor- 
bates on a surface from the diffuse LEED (DLEED) pat- 
terns in which the electron intensity between the Bragg 
spots is capable of a holographic interpretation^. Recent 
work has shown that essentially the same numerical al- 
gorithm is able to form an image-.of the local geometry of 
Si on a SiC(lll)-p(3x3) surfaced. This work was impor- 
tant in that it showed that the holographic method can 
be applied to the case of an ordered array of adatoms on a 
surface by operating on the intensities of superstructure 
(or fractional-order) Bragg spots, and in that it actually 
was the catalyst to the solution of this rather complex 
structure that was most likely beyond the scope of the 
current techniques of conventional .LEED analysis. 

As has been pointed out earlierEJTJ, a set of super- 
structure Bragg spots from an ordered array of adsor- 
bates may be thought of as sampling the corresponding 
diffuse LEED intensities at the positions of these spots. 
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In the case of a large adsorbate unit cell like the (3x3) 
one in the above cited case, the diffuse intensities may 
be thought as being sampled over a rather dense grid in 
the reciprocal-space planes parallel to the surface. How- 
ever, within the class of ordered adsorbate superstruc- 
tures, those of smaller unit cells are much more frequently 
observed. One such set of superstructures which still 
can present difficulties for present state-of-the-art con- 
ventional LEED methods are the p(2x2) structures. It 
is the aim of the present paper to show that even this 
class of adsorbate superstructures may also be brought 
within the purview of holographic LEED methods. 



II. THE 'CORRECT' ALGORITHM FOR 
DIFFUSE LEED 

Holographic methods have already been employed suc- 
cessfully to reconstruct the local geometxx of disordered 
atomic adsorbates on a crystal surfaceEdnl3. In the corre- 
sponding diffuse LEED patterns, diffraction intensity can 
be measured in the areas between the Bragg spots. This 
intensity arises from elastic scattering from a surface ex- 
hibiting structural elements lacking lateral translational 
symmetry. If one considers elastic scattering from a sur- 
face containing a single adsorbate atom, it is obvious that 
the DLEED intensity H a it gives rise to can arise from 
just those scattering paths that include a scattering from 
the adsorbate atom since this is the only one that breaks 
the lateral periodicity of the surface. The diffuse inten- 
sity from a disordered layer of such adsorbates iu. the 
same local sites of a flat substrate may be writtcr£3: 



H(k) =H tt (k)S(ki|) 



(2.1) 



where k = (k|| , k±) is the wavevector of the detected elec- 
tron, k|| its component parallel to, k± that perpendicular 
to the surface, and S is a lattice factor that quantifies 
the degree of long-range order amongst the adsorbates. 
In the case of perfect lattice-gas disorder, which might 
be expected at very low adsorbate coverages, S becomes 
constant and equal to N, the total number of adsorbates 
illuminated by the electron beam, except in the parts of 
the diffraction pattern occupied by the substrate Bragg 
spots (which are excluded in practice), where it is equal 
to TV 2 . Thus the accessible part of the DLEED pattern 
is just a more intense version of H a - 

The holographic view of DLEEDB focuses on the fact 
that the obligatory scattering at an adsorbate leads to 
a natural separation of all scattering paths: electrons 
whose final scattering is by an adsorbate form the refer- 
ence wave -R(k), while those scattered subsequently by 
substrate atoms before reaching the detector provide the 
object wave O(k). On this picture of the adsorbate as a 
microscopic beam-splitter, H a can be written: 



For a single DLEED pattern, it was suggested that 
the holographic reconstruction algorithm derived by 
Bartontj for photoelectron holography could be used to 
produce a 3D image of the local atomic environment of 
the adsorbate. Consequently, the numerical inversion al- 
gorithm can be formulated as a phased 2D Fourier trans- 
form of the DLEED data over k|| . 

Due to the fact that DLEED patterns of different elec- 
tron energies can be measured easily, the idea emerged 
that information from several such patterns could prof- 
itably be combined in reconstructing the-iMagea. Indeed, 
multi-energy reconstruction algorithms 2 !!^ led to a con- 
siderable improvement in the images that could at that 
time be called reliable for the first time. Apart from 
suppressing the unwanted holographic twin image, the 
additional stationary-phase condition arising from the in- 
tegral over energies was designed to single out the contri- 
butions due to the kinematic object wave. A benefit from 
this is that corrections for the remaining anisotropics of 
the reference and object waves may then be performed 
by assuming rather simple forms for these quantitiescj. 
Subsequently, many variants of multi-energy reconstruc- 
tion algorithms have been suggested! 2 !! 

The most successful of such holographic DLEED in- 
version algorithms to date has been the Compensated 
Object- and Referencewave Reconstruction by <m Energy- 
dependent Cartesian Transform (CORRECT)ci. Unlike 
previous multi-energy algorithms that perform the 3D 
holographic reconstruction integral in a polar coordinate 
system (angle and energy), this, transform operates on a 
Cartesian data input (k||, fcjjcll. Considering that a su- 
perstructure spot is characterized by a constant k|| under 
changes of electron energy, this algorithm has the advan- 
tage of being most suitable for an extension to diffraction 
data from ordered superstructures. With this algorithm, 
the reconstructed real space distribution around the ad- 
sorbate |i?(r)| 2 (where r = (rii,z) is a position vector 
relative to an origin at an adsorbate with components 
r* |l parallel and z perpendicular to the surface) can be 
calculated via: 



B{v) = 



[ K(k l{ ,kr,r)xQt\\,k ± )e-( ikr - k ^dk ± 

Jk ± 

(2.3) 



e ik ii- r »d 2 k| 



Apart from the obvious fact that this involves a 3D in- 
tegral over reciprocal space, two noteworthy features are: 
first, it operates not directly on the measured intensities 
but rather on a contrast-enhancing and normalizing 
function 



with 



H a (k) = \R(k)+Q(k)\ 



(2.2) 



H av (kn) 



g(k||,fc ± )-g n „(k||; 

= J H(\n h k ± )dk ± 

Jdkj_ 



(2.4) 



(2.5) 
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It has been shown theoreticallyB that the use of such a 
X-function helps to partially remove the self-interference 
terms |i?(k)| 2 and 0(k)| 2 in the DLEED intensity in Eq. 



(2.2), which give rise to high intensities around the ori- 
gin of the reconstructed real space distribution. More 
importantly, \ has been designed to remove effects due 
to partial ordering amongst the adsorbates. Since for 
fiat substrates (i.e. not those exhibiting steps) the lat- 
tice factor S in Eq. (2.1) is an exclusive function of k||, 
it is unaffected by the integrals over k± in the defini- 
tion of x an d hence cancels out. Deviations from perfect 
lattice-gas disorder at higher coverages, which give rise 
to modulations of S, can thus be efficiently suppressed!!!. 

The anisotropic scattering of low energy electrons by 
the adsorbate beam splitter makes the local environ- 
ment of the adsorbate show up only within the forward- 
scattering cone of the incident wave when images are 
reconstructed by previous inversion algorithmsEl Full 
3D images are only formed with such algorithms by tak- 
ing sets of DLEED data of different electron energies for 
each of a number, (at least two) of different directions of 
incident electronsfcauJ. Of course, this causes consider- 
able experimental effort and so it would be much prefer- 
able to use a set of DLEED patterns from just the most 
reliably-measured normal incidence data which allow to 
improve the data quality by off-line averaging according 
to the underlying surface symmetry. This becomes pos- 
sible with the use of the CORRECT algorithm, due to 
its compensation for the anisotropy of the reference wave 
by estimating its value at the position of a scatterer, by 
inclusion of the kernel 



K(k h k ± ;r) 



(2.6) 



under the integral in (2.3). Here / Q (k; • f) is the atomic 
scattering factor of the adsorbate, k^ the direction of elec- 
tron incidence, and C the so called kernel constant (which 
we take to be real), and which represents a zeroth-order 
approximation to the backscattering by the substrate 
prior to scattering by the adsorbate. An order of mag- 
nitude estimate can be made for C from its theoretical 
derivation (leading to CoNi — 0.75A and Ckn% — 2.20A 
for both the systems O/Ni and K/Ni treated in this 
paper)cJ. In general, the effect of C is well defined and 
hence its value can also simply be optimized, such that 
all atoms in the local adsorption geometry show up with 
similar brightnesses. 

This algorithm-has been shown to give reliable images 
using_thj3oreticalc3'E3, as well as experimental DLEED 
dataBEI. 



III. DLEED AND LEED INTENSITIES 

Even though research on surfaces giving rise to diffuse 
LEED intensities can provide important information on 



initial (disordered) stages of adsorption, the great major- 
ity of structures investigated manifest long-range order in 
the plane parallel to the surface (henceforth referred to as 
the lateral plane). Particularly, many systems, where this 
order results in the surface having larger lateral unit cells 
than that of the bulk, put rather heavy demands on the 
standard quantitative LEED analysis. The large number 
of structural parameters to be determined increases expo- 
nentially the number of trial structures to be considered. 
For such structures, the lattice factor S becomes a sum 
of 5-functions with peaks at the reciprocal lattice rods 
due to the period of the new superstructure. In other 
words, destructive interference between the waves origi- 
nating from different adsorbate-substrate clusters extin- 
guishes all diffraction intensity but that concentrating in 
the newly formed, sharp superstructure spots. 

Nevertheless, these superstructure spots have been 
proven to contain the same crystallographic information 
on the local environment of the adsorbate: as long as 
scattering between different adsorbates can be neglected, 
the diffuse intensity from an adsorbate and its local sur- 
roundings corresponding to a particular value of k|| has 
the same energy dependence as the superstructure spot 
intensity from an ordered array of such adsorbates in an 
equivalent local adsorption geometryt3 and of a surface 
reciprocal lattice vector equal to the same value of k|| . 
As pointed out earlienlH, the superstructure spots may be 
thought of as sampling the DLEED intensity distribution 
of the corresponding lattice gas on a finite grid. It was 
also shown that the neglect of intra-adsorbate scatter- 
ing holds even for relatively dense superstructures under 
normal electron incidenceEa. 

With this understanding, a DLEED holographic in- 
version algorithm like CORRECT may in principle be 
applied to superstructure spot intensities - in particular, 
since the latter takes its input diffraction intensities on 
a Cartesian grid in reciprocal space. The only appar- 
ent difference is that the finite sampling breaks down the 
integral over k|| to a discrete sum. However, a numeri- 
cal implementation of the reconstruction algorithm does 
this even in DLEED, so that the reaL difference is the 
strongly reduced data resolution in k||E3. In this respect, 
the larger the new surface period, the more fractional or- 
der spots and hence the more data left for the algorithm 
and the more similar the situation becomes to the former 
DLEED case. Recently, the CORRECT algorithm was 
thus successfully applied to the p(3x3) reconstruction of 
SiC(lll)o. The high density of superstructure spots cor- 
responding to such a large reconstructed surface unit cell, 
allowed the application of the DLEED algorithm without 
any modifications caused by the finite sampling. Differ- 
ences that appear for small ordered superstructures with 
their even further restricted data base will be addressed 
in the next section. 

Note that the application to SiC(III) highlights the 
fact that for ordered superstructures, it becomes irrel- 
evant for holography whether the microscopic beam- 
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splitter is an adsorbate on top of a substrate or be- 
longs intrinsically to the substrate itself as in the case of 
adatoms in surface reconstructions. This increases con- 
siderably the number of systems to which holographic 
LEED may be applied. The main requirements are that 
the beam-splitter atom is the principal cause of the break 
of the bulk lateral periodicity and that there be only 
one such atom per superstructure unit cell (in order to 
prevent intermixing of images centered on several such 
holographic reference- wave sources). If large, inequiva- 
lent adsorbate overlayer domains exist on the surface the 
resulting image would be a simple superposition of the 
local atomic structure around each adsorbate. Since the 
reconstructed image reveals essentially the local bonding 
geometry of the adsorbates, even the existence of more 
than one inequivalent superstructure domain is immate- 
rial, provided the local geometries of the adsorption sites 
are identical. 

As a last point, it is important to mention that for 
the application to conventional LEED data, the strong 
reduction in data resolution goes hand in hand with 
a strong increase in the available energy range and data 
quality. Conventional LEED spot intensities can be mea- 
sured up to much higher energies (~ 400-500eV and 
even higher at low temperatures) than DLEED distribu- 
tions, which suffer from disturbing effects from the much 
brighter substrate Bragg spotstil and greater contribu- 
tions from thermal diffuse scattering at higher energies. 
Since atomic scattering factors show a much weaker en- 
ergy dependence at elevated energies, this additionally 
accesible data range is also particularly favourable for 
any holographic inversion algorithm, which tend to rely 
on rather smoothly varying electron scattering proper- 
ties. Further, the signal-to-noise ratio of the bright dis- 
crete spots is much higher than that of diffuse intensi- 
ties, and contributions from thermal diffuse scattering 
are less important and easier to subtract. In conclusion, 
the data acquisition of discrete LEED diffraction spots is 
thus much more standard and reliable than that of dif- 
fuse LEED distributions. Consequently, the extension of 
holographic inversion algorithms to ordered systems will 
make this technique much more applicable to a wider 
range of systems. 



the coarse sampling must evidently result in perturbing 
effects on the integral transform designed for continuous 
data distributions. Hence, there are two questions to be 
addressed: (i) Will there be restrictions on the general 
validity of the method? (ii) Is there a way to make more 
efficient use of the smaller quantity of remaining data? 
We will treat question (i) first and defer question (ii) to 
the next section. 

When considering restrictions on the general validity, 
the key point to notice is that the primary consequence of 
the prevailing surface pe riod icity is to reduce the contin- 
uous ky integral in Eq. (2.3) to a discrete sum. Since the 
term inside the square brackets in Eq. (2.3) is a slowly 
varying function of r, this two-dimensional integral is es- 
sentially a Fourier transform in k|| . The problem encoun- 
tered here is hence similar to the well known numerical 
problemj-of discrete Fourier transforms. The sampling 
theoremED states in this case, that when a function that 
is not-bandwidth limited is sampled on discrete intervals 
of magnitude A, there will be a frequency range in the 
dual space limited by the Nyquist frequency / c t3: 



[-fc, fc] with f c 



2tt 
2A' 



(4.1) 



All power spectral density of the original function that 
would lie outside of this range, will be spuriously moved 
inside by aliasing, and the Fourier transform will be pe- 
riodic due to aliasing of spectral density outside of f c . 

For the application of this theorem to our present prob- 
lem, consider for simplicity a laterally one-dimensional 
example with a surface reconstruction with a period of 
n times the bulk one. This gives rise to (n — 1) frac- 
tional order diffraction spots per Brillouin zone in the 
LEED pattern. The sampling interval A will thus be 
g/n, where g is magnitude of the reciprocal lattice vector 
corresponding to the lateral bulk periodicity given by 



<J 



2tt 

a p 



(4.2) 



and a p the lateral lattice constant. Substitution in (4.1) 
determines the Nyquist critical frequency to be 



fc = 



(4.3) 



IV. LIMITS OF VALIDITY FOR SMALL 
SUPERSTRUCTURES 

Even though the above-mentioned application to 
SiC(lll)-p(3x3) employed the unmodified CORRECT 
algorithm, it can only be seen as an intermediary step 
towards the real extension of DLEED holography to or- 
dered superstructures. There are very few such large 
unit cell reconstructions with the additional requirement 
of having only one prominent atom to serve as a holo- 
graphic beam-splitter. In general, surface structures will 
show smaller periodicities with even less fractional-order 
spots remaining for the inversion algorithm. In this case, 



Let us think of some hypothetical superstructure peri- 
odicities: n — 1 represents same periodicity of bulk and 
surface, and hence results in no fractional-order spots. 
Our holographic method is not applicable here. On the 
other hand, n > 2 leads to a larger number of fractional- 
order spots and a larger usable data density. Any re- 
strictions on the algorithm's validity will be expected to 
affect such a case little, and one can argue that it suffices 
to derive such limits for the worst case, namely n — 2. 

From (4.3) above, f c = a p for this case. However, so 
far we have neglected the fact that the integer order spots 
have to be excluded from the inversion algorithm since 
they are dominated by pure substrate scattering paths. 
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This will double the effective sampling interval A to g 
and consequently halve the critical frequency to a p /2. On 
the other hand, for real laterally two-dimensional cases, 
the effective sampling interval may also lie between these 
two estimates, depending on the periodicity in the other 
lateral direction. For example, consider the LEED pat- 
tern from a p(2x2) superstructure: there will be rows of 
superstructure Bragg spots in the direction of each recip- 
rocal lattice vector at intervals of g/2 and alternate rows 
that include integer-order spots where the superstructure 
spots are found at intervals of g. Cum grano salis we may 
therefore argue, that the critical frequency for this case 
must lie somewhere in the range 



< fc < a p . 



(4.4) 



The sampling theorem now implies that any real space 
intensity represented by the function |_B(r)| 2 in Eq. ( |2.3| ) 
that would normally appear outside a lateral range of 
±/ c , is moved inside by aliasing. A natural consequence 
would be severe image distortions due to these contribu- 
tions. Fortunately, the algorithm is saved by the heavy 
damping of the low energy electrons inside the crystal, 
which reduces strongly any contribution to the diffrac- 
tion intensity from atoms farther from the adsorbate. 
Also, holographic fringes due to more distant atoms will 
be of higher frequency in reciprocal space that may be 
beyond the available data resolution. Thus we may 
deduce that any such aliased contributions inside the 
critical frequency range ±/ c be small compared to the 
real intensities corresponding to true atomic positions 
(in Fourier transform terms this can be rephrased that 
-ff(k) is bandwidth limited). This reasoning is corrobo- 
rated by the fact, that in all previous investigations us- 
ing DLEED data, where no such aliasing effects are ex- 
pected, the only atomic positions recovered belonged to 
the direct local adsorption geometry and no intensity due 
to furtb 
imagesE 

The more important consequence of the sampling the- 
orem in the present case is the folding out of intensity 
outside of the critical range f c . In other words, aliasing 
restricts the range of lateral validity of the algorithm to 



ing atoms appeared in the reconstructed 



P r 

— < f c < a p 



(4.5) 



for p(2x2) structures. Outside of this range one can- 
not be sure whether reconstructed intensity is due to 
real atoms or solely due to aliased contributions. Note, 
that even in this worst case, the aliasing does not destroy 
the algorithm's applicability: the latters' main task is to 
recover atomic positions of the local adsorption geom- 
etry, or the local adatom environment. In almost any 
imaginable case the lateral distance of such neighbour- 
ing atoms is less than a p , and hence still accesible to 
the holographic inversion algorithm. In particular, for 
the case of a hollow-site adsorption on a fcc(OOl) surface, 



which we discuss later, the lateral distance of the 4 near- 
est first layer substrate atoms is a p /2 and thus well within 
the critical range. Even though the very restricted data 
resolution will thus not lead to a breakdown of our holo- 
graphic scheme, we must keep clearly in mind its limited 
range in the direction parallel to the surface. 



V. DROPPING THE x-FUNCTION 

We now consider the question whether the limited 
quantity of remaining data in conventional LEED should 
be treated differently by a holographic algorithm. After 
establishing the equivalence, for our purposes, of DLEED 
and LEED intensities, it is obvious that the general form 
of the 3D integral transform does not need to be changed, 
since the algorithm still targets the same scattering-path 
phases. On the other hand, as mentioned earlier, the pri- 
mary purpose for the design of the contrast-enhancing 
%-function in Eq. fl2.4| ) has been to suppress effects of 
partial ordering amongst the adsorbates, which leads to 
modulations in the lattice factor S. However, for the case 
of an ideally ordered superstructure on a flat surface, S 
becomes a sum of (5-functions with equal value at each of 
the remaining reciprocal lattice rods. Any filtering effect 
of x is thus no longer needed. 

The second advantage of presenting data to the inver- 
sion algorithm in the form of the x~f unc tion has been 
t he partial removal of the self-interference terms in Eq. 
([2.2|) which would be expected to give rise to just low- 
frequency oscillations in the diffraction pattern. The 
principal argument for this is that oscillations in i?(k) 
are due mainly to the reference-object interference terms 
R*(k)0(k) and R(k)Q*(k) and that hence ff „(k||) as 
defined in Eq. ( |2.5| ) represents a rough approximation 
to the remaining self-interference terms in H(k)c£. It is 
clear, that this crude estimate causes a degradation of 
the input diffraction data in the construction of the x~ 
function. This is of no importance in the DLEED case, 
where an abundance of data tends to average out such 
deleterious effects, and allows predominantly the positive 
consequences of the x _ function to be noticedEJ. 

With the filtering effect not required for superstruc- 
tures, and knowing that, due to the limited quantities 
of available conventional-LEED data, it will be advanta- 
geous to use rather the un-degraded measured intensities 
H (k) as input to the reconstruction algorithm, we seek 
an alternative way to suppress the self-interference con- 
tributions to the images. This can actually be achieved 
with a much simpler procedure: since these terms in the 
diffraction intensity are primarily of low-frequency in the 
diffraction data, they would be expected to cause distor- 
tions Tnairily in the part of the reconstructed image near 
its origiroBEaa. Given that the origin is defined by the 
position of the atomic beam-splitter, no other atoms are 
expected within a hard-sphere radius of that adatom or 
adsorbate. We therefore suggest, that ignoring any image 
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features within such a radius (i.e. setting B(r) to zero 
in that region), is as effective a method of suppressing 
the self-interference terms as by the use of a x-function, 
but without the danger of degradation of the useful holo- 
graphic fringes in the diffraction data. 



VI. TEST CASE OF O/Ni(001) p(2x2) 

With these modifications, our algorithm should be 
applicable to LEED data from smaller ordered super- 
structures. Such data are characterized by a lower 
reciprocal-space density of fractional-order spots, and 
there may indeed be a lower limit to this density for 
effective holographic reconstruction. Thus in the or- 
der of surface structures, as characterized by overlayer 
sizes, holographic LEED promises to occupy a niche es- 
sentially complementary to that of conventional quan- 
titative LEED methods: the smaller surface unit cells 
possess a limited number of free structural parameters, 
and can be almost routinely solved by standard trial- 
and-error proceduresQ, while the larger superstructures 
are more problematical due to the exponential scaling 
with the number of parameters, as discussed earlier. 
In contrast, the larger density of fractional-order Bragg 
spots from the larger overlayer unit cells makes their 
LEED pattterns more favorable for an holographic anal- 
ysis, as our earlier work on the SiC(lll) (3x3) has al- 
ready demonstrated^. In the present paper we show that 
even superstructures as small as p(2x2), which have in 
the past posed considerable difficulties for conventional 
LEED analysis, may be brought within the purview of 
the holographic method. 

As an initial test system we chose the p(2x2) phase of 
O/Ni(001). Because this is an extensively-investigated 
and well-known system, it is appropriate for testing 
methodologic advances: since Ni is a strong scatterer 
and O a weak one adsorbed very close to the surface, 
it can be argued that this structure represents a rather 
unfavourable scenario for holographic techniques, which 
rely on dominant scattering by the O atom to enhance 
the reference wave and a suppression of intra-Ni and O-Ni 
multiple scattering. A successful recovery of the adsorp- 
tion site in this case would suggest a general applicability 
of our method, especially on e.g. lighter scattering sub- 
strates like Si. 



VII. IMAGE DISPLAY SCHEME 

The reconstructed images in this paper will be dis- 
played in perspective to reveal their full 3D real space 
geom etry. The function |-B(r)| 2 is calculated from Eq. 
(2.3) on a grid of 0.2 A resolution inside a cylinder of 
depth 3. 5 A and a lateral radius of 2.0 A, consistent with 
our estimate of the limits of lateral validity for the unit 



cell parameter a p = 2.49A for Ni. All atoms that spec- 
ify the local adsorption geometry of O/Ni(001)-p(2x2) 
are located within this volume and their positions should 
therefore be retrievable by our method. Small spheres, 
with diameters proportional to the reconstructed inten- 
sity |-B(r)| 2 and whose gray shading varies from white 
for the lower intensities to black for the higher ones, are 
drawn at the grid points. The origin of the coordinate 
system is defined by the adsorbate (beam-splitter), which 
is artificially added to facilitate the understanding. No 
intensity is calculated inside a sphere of l.oA around this 
origin to suppress the self-interference terms as reasoned 
above. 

Even though the commonly used 2D cut-planes 
through a reconstructed image may be mere suited for 
an exact analysis of the atomic positionscS, we believe 
that the overall 3D presentation of the reconstruction re- 
sult not only is more honest (no choice of presentation 
of "particular" cutplanes), but also focuses better on the 
main task of the holographic technique: a direct insight 
into the essential features of a structure, i.e. in this case 
the adsorption site. 



VIII. SIMULATED DATA FROM 
UNRECONSTRUCTED SUBSTRATES 

As a first step we simulated LEED I(V)-data for a 
simplified O/Ni(001)-p(2x2) geometry: O residing in 
the four-fold symmetric hollow site of a bulk-tmncated 
Ni(001) surface at an adsorption height of 0.9AEI. The 
simulation was carried out with the standard van Hove- 
Tong computer codeel using up to 11 phase shifts repre- 
senting atomic scattering properties (this was found to be 
appropriate for the envisioned energy range). The data 
set generated comprises the I(V)-curves of all fractional 
order spots for the energies 100-400 eV. This represents 
a typical and experimentally feasible range for LEED 
investigations. To closely simulate the conditions of a 
typical LEED experiment, the only data used was that 
collectible from within a 50° polar semi-angle, which cor- 
responds to the standard angular range of conventional 
electron detectorsa, and additionally had been shown to 
be the most appropriate in earlier holographic DLEED 
investigations^. 

The reconstructed image shown in Fig. 1 clearly iden- 
tifies the four-fold adsorption site. An unambiguous dis- 
tinction between the 5 atoms of the Ni local adsorption 
geometry and remaining residual "noise" on the image is 
easily possible. The highest such noise at non-atomic lo- 
cations is at 44% of the intensity of the identifie d at oms. 
The constant C in the integral kernel K of Eq. (2J3) has 
been optimized to a value of 0.70A to give approximately 
equal intensity for both 1st and 2nd layer atomstj. This 
value agrees very well with the order-of-magnitude esti- 
mate derived from theory as described in section II. 

The most important property of a direct inversion algo- 
rithm when applied to a priori unknown systems, is its 
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FIG. 1. Recovered local geometry of the O adsorption site 
from LEED data calculated for a model Ni(001)-p(2x2) sur- 
face with an unreconstructed substrate. The electron energy 
range of the data used is 130-400 eV, the kernel constant 
C — 0.70A and the maximum noise level in the image was 
found to be 44% of the maxima denoting the atom positions. 
For details on the display procedure, see Section VI. 

stability. No fundamentally different information, that 
could lead to wrong predictions of structural features, 
should be reconstructed when changing the parameters 
of the algorithm. With most of the data base parameters 
dictated by the experiment, the most influential one left 
to be adjusted is the energy range of the data. Since the 
effect of the finite energy integral in Eq. ( |2.3| ) is only 
to suppress, and not to erase multiple scattering contri- 
butions, it is to be expected that the image quality can 
vary as a function of the energy range used. Due to 
anisotropics in the atomic scattering factors sometimes 
even smaller ranges may lead to improved images. How- 
ever, changing the upper and lower energy boundaries 
should not result in completely different images, on the 
basis of which wrong adsorption sites would be infered. 
We checked on this and found the image to be stable 
when increasing the lower boundary or decreasing the 
upper one by approximately 60eV, i.e. the maximum in- 
tensities were always found at the site of the 5 identified 
atoms. The use of even smaller energy ranges caused 
the collapse of the algorithm as expressed by very noisy 
images with many new intensity peaks. So, no false ad- 
sorption site could have been obtained from the present 
data. 



z 




FIG. 2. Same as Fig. 1, except that the LEED data is 
calculated from a model that had been fitted to experimen- 
tal LEED data in a previous quantitative LEED study. This 
model includes reconstruction of the substrate and fitted ther- 
mal vibration parameters. The electron energy range used 
was 112. 5-362. 5eV, the kernel constant C = 0.95A and maxi- 
mum noise level found to be 37% of the maxima at the atoms. 

The vertical position of all atoms is retrieved to within 
a 0.2A accuracy. As already experienced with DLEED 
dataEj, the lateral position of the 1st layer Ni atoms, that 
appear only due to the inclusion of the integral kernel K , 
has a much greater uncertainty of ~ 0.5A. Also, this lat- 
eral position can shift inwards or outwards (see results in 
the next section) depending on the exact scattering geom- 
etry. This effect has to be related to the simplicity of the 
compensating kernel, which at its present stage includes 
only a zeroth order appipximation to the backscattering 
properties of the systemc.3 and is hence very sensitive to 
details of the atomic scattering factors. However, the un- 
ambiguous recovery of the essential features of a structure 
(like the adsorption site) , is more than sufficient to allow 
a subsequent refinement by the conventional quantitative 
LEED analysis. 

To verify the generality of the result obtained, we sim- 
ulated a second data set, with a different adsorbed chem- 
ical species and adsorption height, but with all other 
parameters unchanged. Consistent with the determina- 
tion of the related K/Ni(001)-c(4x2) structural, we po- 
sitioned K in the four-fold hollow site with an adsorption 
height of 2. 56 A in a hypothetical p(2x2) arrangement. 
Fig. 2 shows that this adsorption geometry is recon- 
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IX. REALISTIC SIMULATION FROM 
RECONSTRUCTED SUBSTRATES 




FIG. 3. Same as Fig. 1, except that the LEED data is 
calculated from a model that had been fitted to experimen- 
tal LEED data in a previous quantitative LEED study. This 
model includes reconstruction of the substrate and fitted ther- 
mal vibration parameters. The electron energy range used 
was If 2.5-362.5eV, the kernel constant C = 0.95A and maxi- 
mum noise level found to be 37% of the maxima at the atoms. 

structed with comparable quality and clarity as in the 
case of O/Ni. The optimized kernel constant C = 2. 35 A 
again corresponds very well to its theoretical estimate. 
The slightly higher noise level of 63% is probably related 
to the K scattering factor, but nevertheless allows an 
unambiguous identification of all 5 atomsEil. The new 
vertical positions due to the enlarged adsorption height 
are again correct within 0.2A and the lateral shift of the 
1st layer atoms is comparable to the previous case. 

The stability encountered with respect to changes of 
the energy range used is even increased with respect to 
O/Ni (variations of up to 90eV are possible), which cor- 
roborates our hypothesis of O /Ni as an unfavourable sys- 
tem for holography: the smaller multiple K-Ni scattering 
due to the increased adsorption height together with the 
stronger reference wave due to the K, reduce the required 
minimum energy range for multiple scattering suppres- 
sion. Consequently, we are led to believe that the two 
results presented here suggest the general applicability 
of the modified CORRECT algorithm. 



A typical concern often raised about results of the type 
presented in the last section is that such simplified sim- 
ulations cannot be compared to real data as obtained 
from any realistic structure. Since the van Hove-Tong 
computer code takes proper account of all multiple scat- 
tering, the main difference between data as in the last 
section and theoretical data as used for actual quantita- 
tive LEED analyses consists of two points: thermal vibra- 
tions and structural deviations from bulk-like truncation 
like e.g. layer relaxations or buckling. It is then argued 
that both effects might possibly influence the holographic 
algorithm. 

In order to address this question, we simulated a data 
set, where all structural and non-structural parameters 
involved in the dynamical LEED calculation were taken 
exactly as they had been obtained in the previous-quan- 
titative analysis of the O/Ni(001)-p(2x2) systenM. Co- 
incidentally, our previous test system proved again to be 
a rather unfavourable one with respect to both of the 
above-mentioned features: not only is the determined vi- 
brational amplitude of the O of 0.3A unusually high, but 
also the Ni substrate experiences relaxations of the first 
two layer spacings (di2 = 1.80A and e? 2 3 = 1.75A), and a 
relatively strong buckling of O.lOA in the second layer. In 
common with the previous simulations, the full data-set 
comprised again of all fractional-order spots in the energy 
range 100-400 eV and the holographic reconstruction was 
carried out by exactly the same procedure as described 
in the last section. 

Consequently, the result shown in Fig. 3 is particu- 
larly gratifying. The local adsorption geometry is recov- 
ered just as clearly as from the simplified data. Stability 
with respect to changes of the energy range used and the 
accuracy of the recovered atomic positions are found to 
match very well with the ones described in the preceding 
section. Note, however, that the unstable lateral position 
of the 1st layer atoms now appears shifted inward from 
the correct value as had already been pointed out in the 
last section. Interestingly, the overall noise level in the 
image is with 37% even slightly lower for this "tougher" 
more realistic model. Knowing that thermal vibrations 
lead to a reduction of the effective electron coherence 
length, it might actually be argued that this effect helps 
reduce disturbing multiple scattering contributions due 
to atoms further from the adsorbate and hence benefits 
a holographic inversion. 

However, the theoretical objection concerning struc- 
tural deviations from bulk-like positions has to be taken 
more seriously. Any substrate atom whose position de- 
viates from that expected from its ideal lattice position, 
like the buckled 2nd layer Ni atom, could in principle act 
as another conduit for electrons to the fractional-order 
Bragg spots. The standard argument for neglecting these 
contributions in the holographic theory has been that 
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FIG. 4. Experimental LEED pattern 

from a O/Ni(001)-p(2x2) surface for an electron energy of 
100 eV. The inset shows the I(V)-curve of the (| 0) beam. 

atoms that only differ slightly from their lattice positions 
contribute rather weakly to the resulting diffraction in- 
tensity, in comparison to the major rupture of the peri- 
odicity due to the introduction of a new atom, such as an 
adsorbate or adatom. As a further theoretical test, we re- 
peated the simulation for the realistic O/Ni(001)-p(2x2) 
geometry, but this time with the 2nd layer buckling arti- 
ficially doubled to 0.2A. Even for this case, the complete 
adsorption geometry was again retrievable under other- 
wise identical parameters. The only apparent effect of 
the increased buckling amplitude seemed reflected in a 
rise of the overall noise level to 62%. This can be seen as 
a corroboration of the argument that the breaking of the 
bulk periodicity by a deviation of the position of a sub- 
strate atom leads to contributions to the fractional-order 
intensities that may degrade the simple picture of the ad- 
sorbate or adatom as the only path to those fractional- 
order spots. However, the recovery of a good image even 
in this case, clearly suggests that their contribution is 
small and only serves to increase the background noise 
level on the image. 

These results also help explain the success of our earlier 
worku for the reconstructed SiC(lll)-p(3x3) structure. 
In that case the application of our algorithm to fractional- 
order Bragg spot data was able to recover an accurate 
image of the local site of an adatom despite the fact that 
other atoms in the surface unit cell were deviated from 
their positions in an unreconstructed layer. 

X. EXPERIMENTAL DATA 

As a last step we try to invert experimental data of 
the O/Ni(001)-p(2x2) system. After standard crystal 
preparation, the Ni substrate was subjected to an ex- 
posure 2 • 10~ 8 mbar of O for 60sec at 90 K. Subsequent 
annealing at 500 K led to the formation of a sharp p(2x2) 




FIG. 5. Same as Fig. 3, except that the image is com- 
puted from experimental data from a O/Ni(001)-p(2x2) sur- 
face. The energy range used was 90-344 eV, the kernel con- 
stant C = 0.33 A and the maximum noise level was 50% of 
the maxima at the atoms. 

superstructure pattern as evidenced in Fig. 4. The mea- 
surement was performed using the-standard Video-LEED 
system developed in Erlangentrta and included 4-fold 
symmetry averaging according to the expected rotational 
symmetry of the diffraction pattern. The full data set 
comprises the 8 symmetry inequivalent fractional-order 
beams closest to specular reflection in the energy range 
90-344 eV. Good agreement of this enlarged data set was 
found with the three fractional-order spot JlV)-curves 
used in the earlier quantitative LEED studyta. 

Fig. 5 shows the reconstructed local adsorption geom- 
etry. Although the experimental data set is smaller with 
respect to the total energy range and the total number 
of fractional-order beams, we nevertheless find essentially 
the same stability, accuracy and unambiguity of the im- 
age as described for the other data sets. Even on reducing 
the number of beams used to the 5 nearest to the (00) 
beam in the diffraction patterns, the image continued to 
reliably show the 4-fold hollow adsorption site with an 
only slightly degraded overall image quality. 

This final result shows, that no particularly large data 
set needs to be measured to ensure a proper working 
of the reconstruction algorithm. Rather, it is the same 
I(V)-curves that are measured for a quantitative LEED 
analysis, that also provide the input to the holographic 
inversion. Ultimately, it is only this similarity of required 
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data that can make holographic LEED a practically use- 
ful complement to the established quantitative LEED 
analysis at the present time. 



XI. CONCLUSIONS 

We have shown in this paper that a fully three- 
dimensional image of the adsorption site of atoms form- 
ing a p(2x2) overlayer on a metal surface may be re- 
covered from a data set that might typically be mea- 
sured in the usual practice of conventional LEED stud- 
ies, i.e. one from just the most-reliably measured nor- 
mal incidence LEED data. We have demonstrated this 
first for a simulated data set from models of O/Ni(001)- 
p(2x2) and K/Ni(001)-p(2x2) surfaces with unrecon- 
structed substrates. The usual picture of holographic 
LEED from adsorbates on surfaces assumes that the in- 
tensities of the fractional-order spots arise only from scat- 
tering paths that include a scattering at the adsorbate. 
This is not strictly true in the presence of substrate re- 
constructions which also cause deviations from the peri- 
odicity of the bulk lattice. O/Ni(001)-p(2x2) is a par- 
ticularly well-known system with an adsorbate-induced 
substrate reconstruction. We show that realistic LEED 
I(V)-curves calculated from the accepted model of this 
substrate reconstruction also give a clear and unambigu- 
ous holographic image of the local adsorption site of the 
O atom. A similar image is also reconstructed from ex- 
perimental LEED data from this surface. 
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